First-principles modeling of strain 
in perovskite ferroelectric thin films 



OSWALDO DIEGUEZ 
Department of Materials Science and Engineering, 
Massachusetts Institute of Technology, 
77 Massachusetts Avenue, Cambridge MA 02139, USA 

DAVID VANDERBILT 
Department of Physics and Astronomy, 
Rutgers University, 
136 Frelinghuysen Rd, Piscataway NJ 08854, USA 

February 3, 2008 



Abstract 

We review the role that first-principles calculations have played in understanding the efi'ects of 
substrate-imposed misfit strain on epitaxially grown perovskite ferroelectric films. We do so by ana- 
lyzing the case of BaTiOa , complementing our previous publications on this subject with unpublished 
data to help explain in detail how these calculations are done. We also review similar studies in the 
literature for other perovskite ferroelectric-film materials. 



1 Introduction 

Ferroelectrics materials are of technological importance due to their ability to sustain a macroscopic 
polarization that can be switched by the action of an electric field. In particular, the demand for smaller 
industrial components has drawn attention in the last few years to ferroelectrics in film form. These 
can sometimes be grown epitaxially in a coherent fashion with thicknesses beyond tens of nanometers. 
Part of what makes them interesting is that the properties of the thin films can differ significantly from 
those of the corresponding bulk form. Several factors contribute to these differences, e.g., the thickness 
of the film, the electrical boundary conditions imposed on it, the orientation of the substrate upon which 
the film grows, the presence of defects or stoichiometric variations, and the misfit strain imposed by the 
substrate on the film. The reader interested in the science of ferroelectric films can consult the recent 
reviews [Tl[2l[3l|4l[5l[6l[7l[8], especially the one by Rabe [6] that specifically addresses the theory of 
epitaxial strain effects. 

In this paper we will be concerned with the role of the misfit strain in the structural behavior of 
ferroelectric thin films. In doing so we will restrict ourselves to perovskites, one of the most representative 
and studied families of ferroelectric materials. We will review how first-principles theories can be used 
to isolate the effect of epitaxial strain from the other film-related effects mentioned above. Section 
[2] introduces the structural phase diagrams that have proven to be very useful in understanding the 
properties of ferroelectric films. Section [3] explains in detail how first-principles calculations can be used 
to obtain those diagrams, focusing on the case of BaTiOs as a prototypical example of a perovskite 
ferroelectric. This section is based in two of our previous publications [9l[T0], complemented by expanded 
results and detailed explanations that have not previously been published. In Sec. 2] we review similar 
studies that have been carried out for ferroelectric films composed of materials other than BaTiOs. 
Finally, we summarize briefly in Sec. [S) 
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Figure 1: Phase diagrams of epitaxial BaTiOa as predicted by the theory of Pertsev el aL[TT]: (a) Using 
the parameters quoted in Ref . . (b) Using the parameters quoted in Ref. ,12j . First- and second-order 
phase transitions are represented by thick and thin hues, respectively. (Reprinted from Ref. 1^). 



2 Pertsev diagrams 

In a seminal paper, Pertsev, Zembilgotov and Tagantsev [111 introduced the concept of mapping the 
equilibrium structure of a ferroelectric perovskite material as a function of temperature and misfit strain, 
thus producing a "Pertsev phase diagram" of the observable epitaxial phases. The effect of epitaxial strain 
is isolated from other aspects of thin-film geometry by computing the structure of the hulk material with 
homogeneous strain tensor constrained to match a substrate having square surface symmetry and a 
given in-plane lattice constant. In addition, short-circuit electrical boundary conditions are imposed, 
equivalent to ideal electrodes above and beneath the film |11| . Examples of such diagrams are shown in 
Fig. HI the labeling of the phases therein is clarified in Table [TJ Given the recognized importance of strain 
in determining the properties of thin-film ferroelectrics, Pertsev diagrams have proven to be of enormous 
value to experimentalists seeking to interpret the behaviour of epitaxial thin films and heterostructures. 

In Ref. llj, the mapping was carried out using a phcnomenological Landau-Devonshire model taken 
from the literature. Within this kind of theory, a tlicrmodynamical potential describing the behaviour of 
the material is expanded in the relevant degrees of freedom. In the case of Ref. [TT] , the basic degrees of 
freedom were stress and polarization, with the temperature being incorporated via a linear temperature 
dependence of certain expansion parameters. For fixed misfit strain and temperature, a minimization over 
the three components of the polarization was carried out. Depending on the direction the polarization 
points at after the optimization, the film is in a distinctive phase labeled as indicated in Table [1] 

This approach should give excellent results in the regime of temperature and strain in which the 
model parameters were fitted. However, it will generally be less accurate when extrapolated to other 
regimes. As an illustration of this, we pointed out in Ref. 9J that different sets of parameters can lead to 
significantly different diagrams even when exactly the same model is used, as shown in Fig. [TJ Of course, 
such an empirical approach is also limited to materials for which all the needed experimental information 
is available. 



Table 1: Summary of possible epitaxial BaTiOa phases. In-plane cell vectors are fixed at ai = ai, 
Sli = ay. Columns list, respectively: phase label; space group; out-of-plane lattice vector; number of free 
internal displacement coordinates; and form of the polarization vector. 
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3 First-principles calculations: the example of BaTiOs films 



In this section we describe in detail how to apply a first-principles theory to the study of misfit strain 
effects in single-domain perovskite-oxide thin films grown on a substrate with square-lattice symmetry. 
In particular, we focus on the prototypical example of barium titanate. In Sec. 13.11 we show how this is 
done at a microscopic level, considering the physics of the individual electrons and ions in the film; this 
is what we call full first-principles calculations. In Sec. 13.21 we show how to apply a Landau-Devonshire 
theory in which first-principles calculations are used to find the coefficients in a Taylor expansion of the 
thermodynamical potential that describes the behaviour of the film; this is what we call first-principles- 
based calculations. In both cases, the calculations are done for zero temperature. Finally, in Sec. 13.31 
we show the results of applying a refined first-principles-based theory to take into account the effects of 
finite temperature. 

3.1 Full first-principles calculations 

3.1.1 Theoretical details 

First-principles (i.e., ab initio) theories make use of the fundamental laws that govern the behaviour of 
electrons and nucleii to derive the macroscopical properties of a material. Among these, density-functional 
theory (DFT) [T3] as implemented in the Kohn-Sham approach [T3] is the most widely used to deal with 
systems such as those described here, since relatively large systems with many electrons can be treated 
accurately with a modest expenditure of computer time. Nowadays, several sophisticated computer codes 
capable of running such DFT calculations are available as open-source or inexpensive packages, including, 
e.g., ABINIT[l5j, the Quantum-ESPRESSOp] suite of programs, Siesta[17], and VASP[l8]. 

The first-principles DFT calculations described in this subsection were carried out using the VASP[T8] 
software package. The electron-ion interaction was described by the projector augmented-wave method 
|19j : semicore electrons are included in the case of Ba (5s^5p^6s^) and Ti (3s^3p^4s^3(i^). The calculations 
employed the Ceperley- Alder [20] form of the local-density approximation (LDA) exchange-correlation 
functional; other flavours of this functional, like the generalized gradient aproximation, do not lead to 
substantial improvements in the case of BaTiOa [21]. We use a plane- wave basis set with a 700 eV 
kinetic-energy cutoff and a 6 x 6 x 6 Monkhorst-Pack [22] sampling of the Brillouin zone. 

Readers wishing to become more familiar with the terminology used in the previous two paragraphs, 
or to know more about the conceptual framework behind DFT, can consult a recent review on the subject 
by Capelle [23]; see also references therein. 

3.1.2 Search for the most stable phases 

Given the constituents of the unit cell, DFT codes allow one to carry out a search for the structural 
configuration with the lowest energy. We did this for the case of BaTiOs in Ref. [9J where, following 
Pertsev et al. [llj, we worked with a bulk-like periodic system in which the film boundary conditions 
imposted by the misfit strain have been applied. We emphasize that there are no surfaces or interfaces 
in our calculation. Essentially, we take a single unit cell from the deep interior of the film and replicate 
it periodically in all three dimensions; this periodic structure is the one used for our calculations. 

We performed optimizations of distorted configurations of the five-atom unit cell of the cubic per- 
ovskite structure in order to look for the six possible phases mentioned in Table [1] Each starting structure 
had the Ba atoms fixed as determined by the substrate, and we displaced the Ti and O atoms along the 
direction of the polarization vector of each phase. Then, while retaining the symmetry determined by 
these displacements, we relaxed the atomic positions and the out-of-plane cell vector until the value of 
the Hellmann-Feynman forces and zz, yz and zx stress tensor components fell below 0.001 eV/A and 
0.005 eV, respectively. This process allowed us to determine the lowest-energy structure for each of the 
six possible phases. 

Following this recipe for different values of the in-plane strain on the film, we obtained the energy 
curves of Fig. [5] (top panel) . In some cases we started with a structure with a given symmetry, but found 
that the system relaxed to a structure with a higher symmetry; for example, starting with a structure 
having the symmetry of the r phase at misfit strains exceeding ~ 0.6%, we find that it converges to a 
structure of aa symmetry (collapse of the green curve onto the orange curve in the figure). It can be 
seen that the c phase is the most stable structure for large negative (compressive) misfit strains, then the 
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r phase is most stable for an intermediate strain range, and finally the aa phase is most stable at large 
positive (expansive) strains. 

Because the phase transitions from c r ^ aa are both of second order, the curves join each other 
so smoothly that it is difficult to locate the precise phase boundaries from the energy curves alone. The 
boundaries can be located much more precisely by using a stability analysis. At each value of misfit in 
the c phase, for example, we carry out finite-difference calculations of x forces and the xz stress as the 
atomic x coordinates and the xz strain are varied. The lowest eigenvalue of the resulting 6x6 Hessian 
matrix is plotted as the blue curve in the middle panel of Fig. [21 and its zero crossing identifies the critical 
misfit. A similar analysis is used to consider z displacements and shear strains in the a and aa phases 
(red and orange lines respectively). By properly considering the zone-center phonons, the elastic shear, 
and linear cross-coupling between these two kinds of degrees of freedom, this analysis allows us to locate 
the second-order phase boundaries much more precisely than is possible through direct comparison of 
total energies. Since the a phase is never the lowest-energy structure for any misfit strain (see top panel), 
we can ignore it, and we locate the c <-> r and r <-» aa transitions as indicated by the vertical dashed lines 
in the figure. 

With first-principles calculations it is also possible to inspect the individual displacements of the atoms 
in each structure, as shown in the bottom panel of FigO At the phase transition points, displacement 
patterns arise that characterize the phases on both sides of the transition. 

The three panels in Fig. [5] give therefore a consistent picture showing how the structural properties 
of the film vary as a function of the misfit strain imposed on it. For large compressive strains, the lowest 
energy corresponds to the c phase; on the other hand, for large tensile strains the aa phase is favored. At 
a misfit strain of Smax(c) — —6.4 x 10"'^ (a = 3.930 A), there is a second-order transition from the c phase 
to the r phase, with the polarization in the r phase continuously rotating away from the z direction as 
the misfit strain increases. At misfit strain Smin(aa) — 6.5 x 10~'^ (a — 3.981 A), the r phase polarization 
completes its rotation into the xy plane, resulting in another second-order transition to the aa phase. 

3.1.3 Comments on the underlying assumptions 

Here we comment briefly on the effects of the assumptions made in the construction of this first-principles 
Pertsev diagram. In principle, we should consider the possibility of equilibrium structures with larger 
unit cells, particularly those with cell-doubling octahedral rotations (i.e., tilts). Such rotations have 
been shown to be important in SrTiOa, and could condense in BaTiOa under sufficiently large misfit 
strains. As an example, we have checked that the paraelectric phase of the film is stable with respect to 
octahedral rotations about the [001] direction (with M3 symmetry) up to an epitaxial compressive strain 
of —70.9 X 10"'^ (a = 3.675 A), far larger than those likely to be experimentally relevant (see Fig. [3]). 

In addition, while we have studied only the effects of epitaxial strain, other physical effects may also 
be relevant to the structure and properties of thin films, such as atomic rearrangements at the film- 
substrate interface and free surface and the instability to the formation of multiple domain structures 
[24]. Finally, our theory relies on the LDA to compute the exchange and correlation terms in DFT. 
This introduces small systematic errors in the calculation, the most important of which is probably the 
error in the equilibrium lattice constant. However, such errors are well understood and well characterized 
in perovskites, and tend to be similar for different materials of this class, so that there is a tendency 
for cancellation of errors when making relative comparisons of quantities such as misfit strains (see, for 
example, Ref. [25]). 

It should be kept in mind, however, that for both films and superlattices, the assumptions that the 
system is in a single domain and that the epitaxial strain strongly dominates other factors will not be 
valid in all cases. Phase diagrams including multiple-domain states have, for example, been discussed in 
Refs. [26l [27l [24] . Other influences that may be important include surface relaxation and reconstruction, 
atomic and electronic rearrangements at the interface, imperfectly compensated macroscopic electric 
fields, deviations from stoichiometry, and the presence of defects. 

3.2 First-principles-based calculations: a Landau-Devonshire model 

In a Landau-Devonshire approach, a ferroelectric system is described by a thermodynamical potential 
that is expanded as a Taylor series of the relevant degrees of freedom. For example, in the work of Pertsev 
and collaborators [llj . the model is built by starting with the bulk free energy expanded in polarization 
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Figure 2: Top panel: energies of the possible epitaxial BaTiOa phases for different misfit strains, as ob- 
tained from the full first-principles calculations. Center panel: value of the lowest non-trivial eigenvalue of 
the force-constant matrix for the three relevant phases; when it is zero, the second-order phase transitions 
denoted by vertical lines occur. Bottom panel: displacements of the atoms from the cubic perovskite cell 
positions, for the most energetically favorable configuration at a given misfit strain; A3(Ti) labels the 
displacement of the Ti atom in units of the third lattice vector, etc. 
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Figure 3: Value of the restoring force on the oxygen octahedra when they are subjected to a rotation 
about the [001] axis of 0.23° in a pattern of M3 symmetry. A negative value indicates an instability 
leading to a doubled unit cell with oppositely rotated octahedra, but this only happens at a very high 
compressive strain (the equilibrium lattice constant is 3.96 A). 

and stress. The reference state is the paraelectric cubic perovskite phase at the bulk critical temperature 
Tc, and the parameters are fit to reproduce experimental observations of the behavior near the bulk 
ferroelectric transition. For the dependence on epitaxial strain, a Legendre transformation is then made 
to obtain the potential as a function of polarization and misfit strain. Because of the way in which the 
parameters are fit, this Landau-Devonshire potential will give its most accurate results for small misfit 
strains and temperatures near the bulk Tc. 

It is also possible to develop models in the same spirit, but constructed from first-principles calcu- 
lations. Compared with the phenomenological approach, such a first-principles-based theory has the 
advantage that the information needed for the model fit can readily be calculated in a consistent and 
accurate manner within the framework of the chosen first-principles approach. Then, once the fitting has 
been completed, such models can be used to compute detailed properties orders of magnitude faster than 
could be done using the full first-principles calculations. Thus, they can be applied to larger systems 
where the full calculations would be impractical. In this Section, we summarize how one such model 
[lOj can be applied to study BaTiOa, providing additional details that were not included in the previous 
publication in order to clarify how the expansion coefficients were obtained. 

3.2.1 Formalism 

The starting point of this analysis is the parameterized total-energy expression presented by King-Smith 
and Vanderbilt in Ref. j25j . This is a Taylor expansion around the cubic perovskite structure in terms 
of the six independent components rji of the strain tensor (i is a Voigt index, i =1-6) and the three 
Cartesian components Ua (a = x, y, z) describing the amplitude of the soft mode defined by the pattern 
of eigen-displacements associated with the smallest eigenvalue of the (zone-center) force-constant matrix. 
The energy (per unit cell) of the bulk material can be expressed as a sum 

E = E-'--{{j^,}) + E^°'\{u^}) + i?'"*({r,a, {u^}), (1) 

of a term E'^'^^^ arising purely from strain, a term E^°^'^ arising purely from the soft-mode amplitude, 
and a term E™'^ describing the interactions between these two kinds of degrees of freedom. The zero of 
energy is taken to correspond to the cubic reference structure. For crystals with cubic symmetry, the 
strain energy is given, to second order in strain, by 

E^'^'^im}) = \BiM + ril + nl) + B^2{mm + mvz + mm) + \Bii{rf^ + Ti + r,l), (2) 

where i?i2, and i?44 are related to the elastic constants of the crystal by factors of the cell volume. 
The soft-mode energy given in Ref. [25J contains terms up to fourth-order in the soft-mode amplitude, 

i^'°^*({wa}) = nu^ + au^ + liulul + ulul -t- ulul), (3) 
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where v? = + Uy + u^, k is twice the soft-mode eigenvalue, and a and 7 are the two independent 
symmetry-ahowed fourth-order coefficients describing the cubic anisotropy. Finally, the interaction be- 
tween the strains and the soft-mode amplitude is given to leading order by 

+Biyz{'niUyUz H- r]5UzUx + Tjeu^Uy), (4) 

where Bixx, Biyy, and B4yz are the phonon-strain interaction coefhcients. All the coefhcients in these 
three parts of the total-energy expression can be obtained from first-principles calculations on a series 
of distorted structures as described in Refs. [25l[T0]. In Fig. [4] we present previously unpublished results 
of our first-principles calculations that were used as the basis for fitting these coefficients for the case of 
BaTiOg. 

Once we have an expression for the energy of bulk BaTiOa, we can develop a potential appropriate 
for describing films. In the case of coherent epitaxy, where strain elements rji, 772 and rj^ are constrained 
while the others are free, this expression is 

G = {Ar,r-,f + Af^^fja + A^^a^) + {Br-jfj + B^a + B)uly + {CrjV + C^cr ^ 

+Duly + Euj + Fulyul + Huly sin^ 9 cos^ 9. (5) 

where cr = a^^ is a possible external uniaxial stress applied perpendicular to the film, u^y — {u^ + 
and 9 — tan^^ (ux/uy), so that 

Ux = Uxy cos 9, (6) 
Uy — UxySm9. (7) 

The detailed expressions relating the coefficients appearing in Eq. ([5]) back to the KSV coefficients in 
Eqs. dm]) are given in Ref. [TO] . 

Once we have determined the set of coefficients in Eq. ([5]), we can predict the phase diagram as a 
function of misfit strain fj and the normal external stress a by minimizing G to find the values of the 
ground-state soft-mode amplitude components. For a fourth-order theory like the present KSV expression, 
the entire optimization process can be done analytically, since it is possible to compute first and second 
derivatives of G and to do a stability analysis of the various possible phases. 

3.2.2 Testing the model: zero perpendicular external stress 

We first consider the usual case in which the external perpendicular stress a vanishes. Figure [5] (top) 
shows the energy curves of the various phases as predicted by the first-principles-based KSV theory (right 
panel) compared with the full DFT results (left panel) . The agreement between the two sets of results is 
very good, with the small differences arising from two sources. First, the first-principles calculations in 
Ref. [5] (and in Sec. 13.11 were performed using the projector-augmented wave method, [TO] while the first- 
principles calculations used to obtain the KSV coefficients referred here were performed using ultrasoft 
pseudopotentials HH] . Second, there are the intrinsic errors associated with the use of a Taylor expansion 
that has been truncated as described in the previous section; these errors are expected to grow as the 
strain and soft mode amplitudes increase. The bottom panels in Fig. [5] show the displacements of the 
atoms from their centrosymmetric perovskite positions as the strain is varied. Again, the agreement 
between the KSV results (right) and the full DFT results (left) is very good. In particular, the square- 
root behavior predicted by the KSV theory is exhibited by the more exact DFT calculations. 

3.2.3 Other results 

In Ref. [lOj we studied the effect of external stress a applied perpendicular to the film in addition to the 
effect of misfit strain. We did this for BaTiOa and other seven perovskites, constructing a phase diagram 
for each in the two-dimensional space of normal stress and misfit strain. We found that the stress-strain 
phase diagrams obtained for all of the perovskites show a universal topology with straight-line phase 
boundaries meeting at a single crossing point. A detailed interpretation of the those results was given in 
Ref. [To], together with an analysis of the variables that drive each kind of behavior seen for each film. 
In the same article we also presented the computed polarization of the films and related these results to 
the concept of polarization matching. 
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Figure 4: Results of first-principles calculations used to fit the coefficients of the first-principles-based 
model for bulk BaTiOs as described in this section. In each panel, each symbol corresponds to a full first 
principles calculation, while each curve is a low-order polynomial fit to the data. The ground-state lattice 
constant and bulk modulus are obtained from the data in (a), through the use of the Birch equation. Bn 
and i?44 are related to the curvature of the third-order polynomials of graphs (b) and (c); B12 is obtained 
from the relation B = Bn -\- 2Bi2- The coefficient k is half the lowest eigenvalue of the force constant 
matrix, obtained using finite differences and the ability of the first-principles programs to compute forces 
on the atoms. The coefficients a and 7 are obtained from the double- well sixth-order polynomial in (d). 
Finally, Bi^x, B^yy, and B^y^ are proportional to the first derivative at zero strain of the third-order 
polynomials of (e), (f), and (g). 
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Figure 5: Top panel: comparison of BaTiOs energy curves for the six epitaxial phases studied, as obtained 
from direct DFT calculations [9] (left), and from the KSV theory (right). Energies are relative to the 
paraelectric structure at zero misfit strain. The lines in the left panel and the symbols in the right panel 
are provided as guides to the eye. Bottom panel: Comparison of BaTiOs atomic displacements for the 
most stable phase at each given value of strain, as obtained from direct DFT calculations [9] (left), and 
from our KSV model (right). Az{Ti) indicates the displacement of the Ti atom along the z direction, etc. 
Symmetry implies that Ay(Ti) = A,(Ti), A,(02) = Ay(Oi), A„(02) = A,(Oi), A,(02) = A,(Oi), 
and Ay(03) = A2;(03). The lines in the left panel are guides to the eye. 



9 



In-plane Lattice Constant (A) 




-15 -10 -5 5 10 15 
Misfit Strain (10"^) 



Figure 6: Phase diagram of epitaxial BaTiOs obtained using the effective Hamiltonian of Zhong, Van- 
derbilt and Rabe ^29j i30j . All transitions are second-order. 

3.3 Finite-temperature calculations 

The calculations presented so far were all done at zero temperature. In Ref. [9j, we extended our study of 
epitaxial BaTiOs to finite temperatures by using the effective Hamiltonian approach of Zhong, Vanderbilt, 
and Rabe.p9l[30) This method follows the spirit of the first-principles-based calculations described above, 
but now the model is expanded to allow for different distortions in each individual unit cell. The relevant 
degrees of freedom are taken to consist of a ferroelectric local-mode vector in each cell; a displacement 
vector in each cell; and the global homogeneous strain variables. The ferroelectric local-mode and the 
displacement variables describe local (in general inhomogeneous) polarizations and strains, respectively. 
The number of parameters needed to describe such a model is larger, since one needs to include new 
coefficients that quantify the interactions between local-mode variables on neighboring sites. However, all 
of the coefficients in the expanded model are again obtained from fitting to an appropriately constructed 
database of first-principles results. Finally, the effective-Hamiltonian model obtained in this way can 
be used as the basis for finite-temperature Monte Carlo (MC) simulations that can be used to map out 
finite-temperature phase diagrams. 

It is straightforward to impose the constraint of fixed in-plane strain by fixing three of the six elements 
of the homogeneous strain tensor during the MC simulations. For each value of in-plane strain, MC 
thermal averages are obtained for the unconstrained components of the homogeneous strain and the 
average polarization, and phase transitions are identified by monitoring the symmetries of these quantities. 
The resulting Pertsev diagram is shown in Fig. [Gj it is reproduced from Ref. , where additional technical 
details were given. 

The Pertsev diagrams of Figs, [l^a), [Ijb), and [6] share the same topology above and just below Tq: 
p at high temperature, c at large compressive misfit, aa at large tensile misfit, and a four-phase point 
connecting these phases with the r phase at Tq- At lower temperature, there is a drastic difference 
between Figs.[lja) and[ljb), with our theory supporting the latter. While our theory underestimates the 
temperature of the 4-phase crossing point in Fig. [6] by about 100 °C, this is the price we pay for insisting 
on a first-principles approach; indeed, this effective Hamiltonian underestimates the temperature of the 
bulk cubic-to-tetragonal transition by about the same amount. 

4 Studies of other perovskite films 

In this section we review work on the effect of misfit strain on perovskite films other than BaTiOs using 
first-principles methodologies similar to the ones described in the previous section. 
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4.1 KNbOa 



KNbOs is a perovskite whose bulk form is very similar to that of BaTiOa. It displays the same sequence 
of transitions with decreasing temperature: cubic to tetragonal to orthorhombic to rhombohedral, with 
five-atom unit cells. 

Using our zero-temperature first-principles-base approach, we found (lOj that KNbOs and BaTiOs 
also behave in a very similar way in terms of their behavior as a function of epitaxial strain in the context 
of film geometries. KNbOa shows the same c ^ r ^ aa sequence of transitions as does BaTiOa when 
varying the misfit strain at zero temperature in the absence of external perpendicular stress. The first 
transition occurs at a misfit strain of —4.8 x 10"'^, and the second occurs at 5.5 x 10~^. 

4.2 PbTiOs 

As in the cases of BaTiOa and KNbOa, bulk PbTiOa has a ferroelectric ground state with a five-atom 
unit cell. In this case, however, there is only one transition, from cubic to tetragonal, as the temperature 
is reduced. 

A very similar approach to the full first-principles one presented in Sec. 13. II was used by Bungaro and 
Rabe to study the infiuence of misfit strain on PbTiOa films [31]. As in the case of BaTiOa, they find 
a c — > r — s- aa sequence of phase transitions, although in this case the range of stability of the r phase 
is narrower than for BaTiOa (from -t-6.6 x 10"'^ to -1-14.8 x 10~^), and occurs only at expansive strains. 
In the same paper, these authors analyze a superlattice with alternating layers of PbTiOs and PbZrOa, 
finding again the c ^ r ^ aa sequence of transitions, but with a wider region of r behavior than in the 
pure PbTiOs case. 

Our first-principles-based approach of Ref. fTOj also predicts that PbTiOs adopts the c phase for 
compressive strains and the aa phase for tensile strains, with a window of around 1% misfit strain in 
between (from —3.0 x 10^^ to 8.4 x 10~'^). However, we find that the r phase is less energetically favorable 
than a combination of c and aa domains. 

4.3 SrTiOg 

This perovskite is a "quantum paraelectric" in bulk form, meaning that it is only the zero-point motion 
of ions that prevents ferroelectricity from developing. The cubic phase undergoes a non-ferroelectric 
oxygen-tilting (or antiferrodistortive) transition at about 105 K. 

Antons and coworkers [32j carried out a full first-principles investigation along the lines we have 
described in Sec. 13.11 for SrTiOs films with five atoms in the unit cell. They found that in this case the 
sequence of second-order transitions is c —> p ^ aa as the strain goes from compressive to tensile. The 
paraelectric phase exists for misfit strains between —7.5 x 10""^ and -1-5.4 x 10"'^. They also observed a 
strong dependence of the dielectric properties of these films on the misfit strain. A similar behavior was 
observed by us in flO^, although with a narrower window olp phase in between the c and aa phases. 

In their first-principles study, Lin, Huang, and Guo [33] considered larger unit cells for SrTiOs, 
including in this way the possibility of having antiferrodistortive instabilities in the films. As in the 
previous studies, they find that the polarization goes from pointing in the [001] direction for sufficiently 
compressive strains to pointing in the [110] direction for sufficiently tensile strains. The paraelectric 
window extends in their case from —4 x 10~^ to +4 x lO^'^. Their phase diagram is richer due to the 
presence of the new phases that involve rotations of the oxygen octahedra. 

4.4 KTaOg 

As in the case of SrTiOs, ferroelectricity is suppressed in bulk KTaOs by quantum fluctuations. In 
Ref. [34^ Akbarzadeh and coworkers use a first-principles based method to draw the Pertsev diagram for 
KTaOs films. Unlike the theories described above, this work includes electrostatic boundary conditions 
and finite-size effects, albeit in an approximate way. For films that are 28 A thick, these authors find 
that quantum fiuctuations have almost no effect when the electrodes are ideal (short-circuit boundary 
conditions). In this case, the diagram is very similar to that of strong ferroelectrics like BaTiOs. When 
they simulate more realistic electrodes, they find that the corresponding depolarizing fields in the film 
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strongly affect the phase diagram. Under these conditions, the quantum fluctuations might affect the 
phase diagram and eventually suppress the r phase. 

4.5 Other perovskites 

Bulk CaTiOs experiences a phase transition from the cubic phase (with a five-atom unit cell) to an 
orthorhombic phase (with a twenty-atom unit cell) as temperature is reduced. Our calculations TO] for 
CaTiOa films show a similar behavior as for the PbTiOa ones, with the window between c and aa phases 
(misfit strains of —2.3 x 10^^ to 5.35 x 10^^) showing a mixture of c and aa domains. 

The ground state of bulk NaNbOa shows a ferroelectric monoclinic phase with twenty atoms in the 
unit cell. The behavior of the film form [lOj is similar to that of KNbOs, with a c ^ r ^ aa sequence of 
transitions occurring at misfit strains of 5.5 x 10^'^ and 4.1 x 10""^. 

Bulk PbZrOs undergoes a transition from cubic to an antiferroelectric phase that has forty atoms in 
the unit cell. In Ref. [10] we show that the five-atom unit cell film stays in the r phase except for large 
misfit strains that are experimentally difficult to achieve. The behavior of the true ground-state structure 
with epitaxial strain has yet to be fully explored. 

BaZrOs maintains the simple cubic structure in bulk form at all temperatures. According to our 
previous study TOj, it also remains paraelectric in film form for all experimentally relevant misfit strains, 

5 Summary 

In this article we have reviewed the role that first-principles calculations have played in understanding 
the effects of substrate- imposed misfit strain on epitaxially grown perovskite ferroelectric films. To do 
so we have analyzed the case of BaTiOa, complementing our previous publications on this subject with 
unpublished data to help explain in detail how these calculations are done. In particular, we have added 
new explanations and figures to clarify how we carried out a stability analysis to identify precisely the 
critical misfit strains at which phase transitions occur, how we checked for the possible effects of octahedral 
rotations, and how we computed the coefficients of our first-principles based Landau-Devonshire theory. 
We have also reviewed similar studies in the literature for other perovskite ferroelectric-film materials. 

The use of first-principles calculations to understand the effects of misfit strain on the properties of 
films is not restricted to the type of materials discussed here. For example, Fennie and Rabe [3S] have 
suggested that misfit strain can be used as a mechanism to obtain strong coupling between ferroelectric 
and magnetic ordering in EuTiOs. Even more recently, Ishida and Licbsch [3(? have found that, when 
grown epitaxially on SrTiOs, LaTiOa becomes a highly correlated metal, instead of being a Mott insulator 
as found in the bulk form. 

We have shown that important insights can be gained by isolating the effects of misfit strain on the 
properties of ferroelectric films. In the future, it would be desirable to extend the theory by developing 
similar systematic approaches capable of handling other effects that are relevant for films, including 
electric boundary conditions, domain formation, surface and interface effects, and the role of vacancies 
and impurities. 



This work was supported in part by ONR Grant No. N00014-05- 1-0054. 
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